High-fidelity, hyper-accurate, and evolved mutants rewire atomic-level communication in CRISPR-Cas9.

The high-fidelity (HF1), hyper-accurate (Hypa), and evolved (Evo) variants of the CRISPR-associated protein 9 (Cas9) endonuclease are critical tools to mitigate off-target effects in the application of CRISPR-Cas9 technology. The mechanisms by which mutations in recognition subdomain 3 (Rec3) mediate specificity in these variants are poorly understood. Here, solution nuclear magnetic resonance and molecular dynamics simulations establish the structural and dynamic effects of high-specificity mutations in Rec3, and how they propagate the allosteric signal of Cas9. We reveal conserved structural changes and dynamic differences at regions of Rec3 that interface with the RNA:DNA hybrid, transducing chemical signals from Rec3 to the catalytic His-Asn-His (HNH) domain. The variants remodel the communication sourcing from the Rec3 α helix 37, previously shown to sense target DNA complementarity, either directly or allosterically. This mechanism increases communication between the DNA mismatch recognition helix and the HNH active site, shedding light on the structure and dynamics underlying Cas9 specificity and providing insight for future engineering principles.

Substrate-independent activation pathways of the CRISPR-Cas9 HNH nuclease.

A hallmark of tightly regulated high-fidelity enzymes is that they become activated only after encountering cognate substrates, often by an induced-fit mechanism rather than conformational selection. Upon analysis of molecular dynamics trajectories, we recently discovered that the Cas9 HNH domain exists in three conformations: 1) Y836 (which is two residues away from the catalytic D839 and H840 residues) is hydrogen bonded to the D829 backbone amide, 2) Y836 is hydrogen bonded to the backbone amide of D861 (which is one residue away from the third catalytic residue N863), and 3) Y836 is not hydrogen bonded to either residue. Each of the three conformers differs from the active state of HNH. The conversion between the inactive and active states involves a local unfolding-refolding process that displaces the Cα and side chain of the catalytic N863 residue by ∼5 Å and ∼10 Å, respectively. In this study, we report the two largest principal components of coordinate variance of the HNH domain throughout molecular dynamics trajectories to establish the interconversion pathways of these conformations. We show that conformation 2 is an obligate step between conformations 1 and 3, which are not directly interconvertible without conformation 2. The loss of hydrogen bonding of the Y836 side chain in conformation 3 likely plays an essential role in activation during local unfolding-refolding of an α-helix containing the catalytic N863. Three single Lys-to-Ala mutants appear to eliminate this substrate-independent activation pathway of the wild-type HNH nuclease, thereby enhancing the fidelity of HNH cleavage.

High-Fidelity, Hyper-Accurate, and Evolved Mutants Rewire Atomic Level Communication in CRISPR-Cas9.

The Cas9-HF1, HypaCas9, and evoCas9 variants of the Cas9 endonuclease are critical tools to mitigate off-target effects in the application of CRISPR-Cas9 technology. The mechanisms by which mutations in the Rec3 domain mediate specificity in these variants are poorly understood. Here, solution NMR and molecular dynamics simulations establish the structural and dynamic effects of high-specificity mutations in Rec3, and how they propagate the allosteric signal of Cas9. We reveal conserved structural changes and peculiar dynamic differences at regions of Rec3 that interface with the RNA:DNA hybrid, transducing chemical signals from Rec3 to the catalytic HNH domain. The variants remodel the communication sourcing from the Rec3 α-helix 37, previously shown to sense target DNA complementarity, either directly or allosterically. This mechanism increases communication between the DNA mismatch recognition helix and the HNH active site, shedding light on the structure and dynamics underlying Cas9 specificity and providing insight for future engineering principles.

Analysis of coordinated NMR chemical shifts to map allosteric regulatory networks in proteins.

The exquisite sensitivity of the NMR chemical shift to local environment makes it an ideal probe to assess atomic level perturbations in proteins of all sizes and structural compositions. Recent advances in solution and solid-state NMR spectroscopy of biomolecules have leveraged the chemical shift to report on short- and long-range couplings between individual amino acids to establish "networks" of residues that form the basis of allosteric pathways that transmit chemical signals through the protein matrix to induce functional responses. The simple premise that thermodynamically and functionally coupled regions of a protein (i.e. active and allosteric sites) should be reciprocally sensitive to structural or dynamic perturbations has enabled NMR spectroscopy, the premier method for molecular resolution of protein structural fluctuations, to occupy a place at the forefront of investigations into protein allostery. Here, we detail several key methods of NMR chemical shift analysis to extract mechanistic information about long-range chemical signaling in a protein, focusing on practical methodological aspects and the circumstances under which a given approach would be relevant. We also detail some of the experimental considerations that should be made when applying these methods to specific protein systems.

A topological data analytic approach for discovering biophysical signatures in protein dynamics.

Identifying structural differences among proteins can be a non-trivial task. When contrasting ensembles of protein structures obtained from molecular dynamics simulations, biologically-relevant features can be easily overshadowed by spurious fluctuations. Here, we present SINATRA Pro, a computational pipeline designed to robustly identify topological differences between two sets of protein structures. Algorithmically, SINATRA Pro works by first taking in the 3D atomic coordinates for each protein snapshot and summarizing them according to their underlying topology. Statistically significant topological features are then projected back onto a user-selected representative protein structure, thus facilitating the visual identification of biophysical signatures of different protein ensembles. We assess the ability of SINATRA Pro to detect minute conformational changes in five independent protein systems of varying complexities. In all test cases, SINATRA Pro identifies known structural features that have been validated by previous experimental and computational studies, as well as novel features that are also likely to be biologically-relevant according to the literature. These results highlight SINATRA Pro as a promising method for facilitating the non-trivial task of pattern recognition in trajectories resulting from molecular dynamics simulations, with substantially increased resolution.

Redox-dependent structure and dynamics of macrophage migration inhibitory factor reveal sites of latent allostery.

Macrophage migration inhibitory factor (MIF) is a multifunctional immunoregulatory protein that is a key player in the innate immune response. Given its overexpression at sites of inflammation and in diseases marked by increasingly oxidative environments, a comprehensive understanding of how cellular redox conditions impact the structure and function of MIF is necessary. We used NMR spectroscopy and mass spectrometry to investigate biophysical signatures of MIF under varied solution redox conditions. Our results indicate that the MIF structure is modified and becomes increasingly dynamic in an oxidative environment, which may be a means to alter the MIF conformation and functional response in a redox-dependent manner. We identified latent allosteric sites within MIF through mutational analysis of redox-sensitive residues, revealing that a loss of redox-responsive residues attenuates CD74 receptor activation. Leveraging sites of redox sensitivity as targets for structure-based drug design therefore reveals an avenue to modulate MIF function in its "disease state."

Structural Basis for Reduced Dynamics of Three Engineered HNH Endonuclease Lys-to-Ala Mutants for the Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-Associated 9 (CRISPR/Cas9) Enzyme.

Many bacteria possess type-II immunity against invading phages or plasmids known as the clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated 9 (Cas9) system to detect and degrade the foreign DNA sequences. The Cas9 protein has two endonucleases responsible for double-strand breaks (the HNH domain for cleaving the target strand of DNA duplexes and RuvC domain for the nontarget strand, respectively) and a single-guide RNA-binding domain where the RNA and target DNA strands are base-paired. Three engineered single Lys-to-Ala HNH mutants (K810A, K848A, and K855A) exhibit an enhanced substrate specificity for cleavage of the target DNA strand. We report in this study that in the wild-type (wt) enzyme, D835, Y836, and D837 within the Y836-containing loop (comprising E827-D837) adjacent to the catalytic site have uncharacterizable broadened 1H15N nuclear magnetic resonance (NMR) features, whereas remaining residues in the loop have different extents of broadened NMR spectra. We find that this loop in the wt enzyme exhibits three distinct conformations over the duration of the molecular dynamics simulations, whereas the three Lys-to-Ala mutants retain only one conformation. The versatility of multiple alternate conformations of this loop in the wt enzyme could help to recruit noncognate DNA substrates into the HNH active site for cleavage, thereby reducing its substrate specificity relative to the three mutants. Our study provides further experimental and computational evidence that Lys-to-Ala substitutions reduce dynamics of proteins and thus increase their stability.

A Cysteine Variant at an Allosteric Site Alters MIF Dynamics and Biological Function in Homo- and Heterotrimeric Assemblies.

Macrophage migration inhibitory factor (MIF) is an inflammatory protein with various non-overlapping functions. It is not only conserved in mammals, but it is found in parasites, fish, and plants. Human MIF is a homotrimer with an enzymatic cavity between two subunits with Pro1 as a catalytic base, activates the receptors CD74, CXCR2, and CXCR4, has functional interactions in the cytosol, and is reported to be a nuclease. There is a solvent channel down its 3-fold axis with a recently identified gating residue as an allosteric site important for regulating, to different extents, the enzymatic activity and CD74 binding and signaling. In this study we explore the consequence of converting the allosteric residue Tyr99 to cysteine (Y99C) and characterize its crystallographic structure, NMR dynamics, stability, CD74 function, and enzymatic activity. In addition to the homotrimeric variant, we develop strategies for expressing and purifying a heterotrimeric variant consisting of mixed wild type and Y99C for characterization of the allosteric site to provide more insight.

Twisting and swiveling domain motions in Cas9 to recognize target DNA duplexes, make double-strand breaks, and release cleaved duplexes.

The CRISPR-associated protein 9 (Cas9) has been engineered as a precise gene editing tool to make double-strand breaks. CRISPR-associated protein 9 binds the folded guide RNA (gRNA) that serves as a binding scaffold to guide it to the target DNA duplex via a RecA-like strand-displacement mechanism but without ATP binding or hydrolysis. The target search begins with the protospacer adjacent motif or PAM-interacting domain, recognizing it at the major groove of the duplex and melting its downstream duplex where an RNA-DNA heteroduplex is formed at nanomolar affinity. The rate-limiting step is the formation of an R-loop structure where the HNH domain inserts between the target heteroduplex and the displaced non-target DNA strand. Once the R-loop structure is formed, the non-target strand is rapidly cleaved by RuvC and ejected from the active site. This event is immediately followed by cleavage of the target DNA strand by the HNH domain and product release. Within CRISPR-associated protein 9, the HNH domain is inserted into the RuvC domain near the RuvC active site via two linker loops that provide allosteric communication between the two active sites. Due to the high flexibility of these loops and active sites, biophysical techniques have been instrumental in characterizing the dynamics and mechanism of the CRISPR-associated protein 9 nucleases, aiding structural studies in the visualization of the complete active sites and relevant linker structures. Here, we review biochemical, structural, and biophysical studies on the underlying mechanism with emphasis on how CRISPR-associated protein 9 selects the target DNA duplex and rejects non-target sequences.

Principles of target DNA cleavage and the role of Mg2+ in the catalysis of CRISPR-Cas9.

At the core of the CRISPR-Cas9 genome-editing technology, the endonuclease Cas9 introduces site-specific breaks in DNA. However, precise mechanistic information to ameliorating Cas9 function is still missing. Here, multi-microsecond molecular dynamics, free-energy and multiscale simulations are combined with solution NMR and DNA cleavage experiments to resolve the catalytic mechanism of target DNA cleavage. We show that the conformation of an active HNH nuclease is tightly dependent on the catalytic Mg2+, unveiling its cardinal structural role. This activated Mg2+-bound HNH is consistently described through molecular simulations, solution NMR and DNA cleavage assays, revealing also that the protonation state of the catalytic H840 is strongly affected by active site mutations. Finally, ab-initio QM(DFT)/MM simulations and metadynamics establish the catalytic mechanism, showing that the catalysis is activated by H840 and completed by K866, rationalising DNA cleavage experiments. This information is critical to enhance the enzymatic function of CRISPR-Cas9 toward improved genome-editing.

The N-terminus of MIF regulates the dynamic profile of residues involved in CD74 activation.

Macrophage migration inhibitory factor (MIF) is an immunomodulatory protein with a pathogenic activity in various inflammatory disorders, autoimmune diseases, and cancer. The majority of MIF-triggered pathological conditions are associated with activation of the cell surface receptor CD74. In the absence of small molecule antagonists that directly target CD74, MIF variants and MIF-ligand complexes have served as modulators of CD74 activity. These molecules have been reported to have either antagonistic or agonistic effects against the receptor, although the mechanistic parameters that distinguish the two groups are largely unknown. Through molecular dynamics simulations and NMR experiments, we explored the relationship between MIF's catalytically active N-terminus and the surface residues important for the activation of CD74. We found that the two sites are connected via backbone dynamics that are propagated to the CD74 activation surface of MIF, from the ß2 and ß4 strands. Our results also provide mechanistic evidence that explain the functional characteristics of MIF variants, serving as CD74 agonists or antagonists. Such findings are of high importance for understanding the MIF-induced activation of CD74 as well as for the development of highly potent CD74 therapeutics.

1H, 13C, 15 N backbone resonance assignment of the recognition lobe subdomain 3 (Rec3) from Streptococcus pyogenes CRISPR-Cas9.

Rec3 is a subdomain of the recognition (Rec) lobe within CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)-associated protein Cas9 that is involved in nucleic acid binding and is critical to HNH endonuclease activation. Here, we report the backbone resonance assignments of an engineered construct of the Rec3 subdomain from Streptococcus pyogenes Cas9. We also analyze backbone chemical shift data to predict secondary structure and an overall fold that is consistent with that of Rec3 from the full-length S. pyogenes Cas9 protein.

High-Throughput Screening of a Functional Human CXCL12-CXCR4 Signaling Axis in a Genetically Modified S. cerevisiae: Discovery of a Novel Up-Regulator of CXCR4 Activity.

CXCL12 activates CXCR4 and is involved in embryogenesis, hematopoiesis, and angiogenesis. It has pathological roles in HIV-1, WHIM disease, cancer, and autoimmune diseases. An antagonist, AMD3100, is used for the release of CD34+ hematopoietic stem cells from the bone marrow for autologous transplantation for lymphoma or multiple myeloma patients. Adverse effects are tolerated due to its short-term treatment, but AMD3100 is cardiotoxic in clinical studies for HIV-1. In an effort to determine whether Saccharomyces cerevisiae expressing a functional human CXCR4 could be used as a platform for identifying a ligand from a library of less ∼1,000 compounds, a high-throughput screening was developed. We report that 2-carboxyphenyl phosphate (fosfosal) up-regulates CXCR4 activation only in the presence of CXCL12. This is the first identification of a compound that increases CXCR4 activity by any mechanism. We mapped the fosfosal binding site on CXCL12, described its mechanism of action, and studied its chemical components, salicylate and phosphate, to conclude that they synergize to achieve the functional effect.

Regulation of MIF Enzymatic Activity by an Allosteric Site at the Central Solvent Channel.

In proteins with multiple functions, such as macrophage migration inhibitory factor (MIF), the study of its intramolecular dynamic network can offer a unique opportunity to understand how a single protein is able to carry out several nonoverlapping functions. A dynamic mechanism that controls the MIF-induced activation of CD74 was recently discovered. In this study, the regulation of tautomerase activity was explored. The catalytic base Pro1 is found to form dynamic communications with the same allosteric node that regulates CD74 activation. Signal transmission between the allosteric and catalytic sites take place through intramolecular aromatic interactions and a hydrogen bond network that involves residues and water molecules of the MIF solvent channel. Once thought to be a consequence of trimerization, a regulatory function for the solvent channel is now defined. These results provide mechanistic insights into the regulation of catalytic activity and the role of solvent channel water molecules in MIF catalysis.

Allosteric Motions of the CRISPR-Cas9 HNH Nuclease Probed by NMR and Molecular Dynamics.

CRISPR-Cas9 is a widely employed genome-editing tool with functionality reliant on the ability of the Cas9 endonuclease to introduce site-specific breaks in double-stranded DNA. In this system, an intriguing allosteric communication has been suggested to control its DNA cleavage activity through flexibility of the catalytic HNH domain. Here, solution NMR experiments and a novel Gaussian-accelerated molecular dynamics (GaMD) simulation method are used to capture the structural and dynamic determinants of allosteric signaling within the HNH domain. We reveal the existence of a millisecond time scale dynamic pathway that spans HNH from the region interfacing the adjacent RuvC nuclease and propagates up to the DNA recognition lobe in full-length CRISPR-Cas9. These findings reveal a potential route of signal transduction within the CRISPR-Cas9 HNH nuclease, advancing our understanding of the allosteric pathway of activation. Further, considering the role of allosteric signaling in the specificity of CRISPR-Cas9, this work poses the mechanistic basis for novel engineering efforts aimed at improving its genome-editing capability.

NMR and computational methods for molecular resolution of allosteric pathways in enzyme complexes.

Allostery is a ubiquitous biological mechanism in which a distant binding site is coupled to and drastically alters the function of a catalytic site in a protein. Allostery provides a high level of spatial and temporal control of the integrity and activity of biomolecular assembles composed of proteins, nucleic acids, or small molecules. Understanding the physical forces that drive allosteric coupling is critical to harnessing this process for use in bioengineering, de novo protein design, and drug discovery. Current microscopic models of allostery highlight the importance of energetics, structural rearrangements, and conformational fluctuations, and in this review, we discuss the synergistic use of solution NMR spectroscopy and computational methods to probe these phenomena in allosteric systems, particularly protein-nucleic acid complexes. This combination of experimental and theoretical techniques facilitates an unparalleled detection of subtle changes to structural and dynamic equilibria in biomolecules with atomic resolution, and we provide a detailed discussion of specialized NMR experiments as well as the complementary methods that provide valuable insight into allosteric pathways in silico. Lastly, we highlight two case studies to demonstrate the adaptability of this approach to enzymes of varying size and mechanistic complexity.